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ABSTRACT 

We present equivalent width measurements and limits of six diffuse interstellar bands 
(DIBs, A4428, A5705, A5780, A5797, A6284, and A6613) in seven damped Lya absorbers 
(DLAs) over the redshift range 0.091 < z < 0.524, sampling 20.3 < logiV(Hi) < 21.7. 
DIBs were detected in only one of the seven DLAs, that which has the highest reddening 
and metallicity. Based upon the Galactic DIB-A(Hi) relation, the A6284 DIB equivalent 
width upper limits in four of the seven DLAs are a factor of 4-10 times below the A6284 
DIB equivalent widths observed in the Milky Way, but are not inconsistent with those 
present in the Magellanic Clouds. Assuming the Galactic T)IB-E(B — V) relation, we 
determine reddening upper limits for the DLAs in our sample. Based upon the E(B — V) 
limits, the gas-to-dust ratios, N(Hi)/E(B — V), of the four aforementioned DLAs are at 
least ~ 5 times higher than that of the Milky Way ISM. The ratios of two other DLAs 
are at least a factor of a few times higher. The best constraints on reddening derive 
from the upper limits for the A5780 and A6284 DIBs, which yield E(B - V) < 0.08 for 
four of the seven DLAs. Our results suggest that, in DLAs, quantities related to dust, 
such as reddening and metallicity, appear to have a greater impact on DIB strengths 
than does Hi gas abundance; the organic molecules likely responsible for DIBs in DLA 
selected sightlines are underabundant relative to sightlines in the Galaxy of similarly 
high A r (Hi). With regards to the study of astrobiology, this could have implications 
for the abundance of organic molecules in redshifted galaxies. However, since DLAs 
are observed to have low reddening, selection bias likely plays a role in the apparent 
under abundance of DIBs in DLAs. 
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PKS 1127-145, Q1229-020) - techniques: spectroscopic 



1. Introduction 



been studied (Galazutdi nov et al 



2000 



Since th eir discovery bvlHegerl (119221) . several hundred diffuse interstellar bands (DIBs) have 



Jenniskens fc Desertlll994l : iTuairisg et alJl200Cl : IWeselak et al 



2000: iHobbs et al.ll2007l ). and yet no positive identifications of the carriers have been made. The 
DIBs span the visible spectrum between 4000 and 13,000 A. Despite no positive identifications, 
several likely organic molecular candidates have emerged as the sources of the DIBs, including 
polycyclic aromatic hydrocarbons (P AHs), fullerenes, long carbon chains, and polycyclic aromatic 



nitro gen heterocycles (PANHs) (e.g.. IHerbig 



1995 



Snow 



2001 



Cox Sz Spaans 



2006a 



Hudgins et al 



2005). The organic- molecular origin of the DIBs may give them an importance to astrobiology; 
PAHs are now considered an imp ortant early constituent to the inventory of organic compounds 
on Earth (IBada Lazcandl2002l ). Via their infrared emitting vibrational bands, PAHs have bee n 
observed in high redshift dusty ultra-luminous infrared galaxies (ULIRGs, e.g.. lYan et al.ll2005l ). 
Searching for DIBs using the technique of quasar absorption lines provides a different approach 
for charting the presence of possible organics to high redshift. As such, observing DIBs in high 
redshift galaxies may offer an independent method for constraining the environmental conditions in 
early-epoch galaxies governing the abundances of organic molecules, determining the cosmic epoch 
at which organic molecules first formed, and ultimately charting their evolution with redshift. 

Aside from the hundreds of detections within the Galaxy (e.g.lGalazutdinov et al hood I Ijenniskens fc Desertl 



1994! : ITuairisg et al.1 1200(1: IWeselak et al 



the Magellanic Clouds (IWeltv et al 



seven s tarburst ga laxies (jHeckman fc Lehnert 



2006 



2000 



Cox et al 

2oooh 



Hob bs et al. 120071') , DIBs have been detected i n 



2006bl . 120071 ). M31 (ICordiner et al.1 120081 ) 



active galaxy Centaurus A via supernova 



1986A (jRichl Il987l ). spiral galaxy NGC 1448 via supernovae 2001el and 2003hn (|Sollerman et al 
20051). one damped Lyq absorber (DLA) at z = 0.524 toward the quasar (QSO) AO 0235+164 
(jjunkkarinen et al.l 12004k [York et al.ll2006l ). and one z = 0.157 Can selected absorber toward QSO 
J0013-0024 ([Ellison et al.1 12OO8I ) ■ 



There are se veral environmenta l factors, such as Hi column density (IHerbig 



1995 



Welty et al 



20061 ). reddening (IWeltv et al.ll2006l ). and metallicity and ionizing radiation (ICox et al.ll2007l ). that 
are rela ted to DIB strengths. In the Galaxy, DIB absorption strengths correlate strongly with 



iV(Hl) (|Herbigill995l ; IWeltv et al.ll2006l ). However, in the Magellanic Clouds, DIBs are weaker by 
factors of 7-9 (LM C) and ~ 20 (SMC) compared to those observed in the Galaxy with similar iV(Hi) 
(jWeltv et al.ll2006l ). It is not known whether other galaxies in the Local Group and beyond obey the 
Galactic DIB-iV(Hi) relation, or, like the Magellanic Clouds, show departures from this relation. 
Galaxies with high iV(Hi) observed in absorption (i.e., DLAs) that reside at low to intermediate 
redshifts (where the prominent DIBs fall in the optical region) provide excellent astrophysical 
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laboratories with which to investigate this issue. 

In this paper we search for A4428, A5780, A5797, A6284, and A6613 DIB absorption in seven 
low to intermediate redshift DLAs. In § [2J we give a brief summary of each intervening DLA in our 
sample. In § El we discuss the spectroscopic observations and data reduction of the background 
QSOs. In § IH we explain the procedure of our analysis, and the resulting spectra. In § [5j we 
present our results and compare our data to the Galactic DIB-A r (Hi) relation, deduce upper limits 
for the reddening, E(B — V), determine lower limits on the gas-to-dust ratios, and discuss the role 
of metallicity for our sample of DLAs. We conclude in § [6j 



2. DLA Sample 



To potentially maximize our chances of detecting DIBs, and to test the Galactic DIB-iV(Hi) 
relation in redshifted galaxies, we selected DLAs toward background QSOs having the highest 
iV(Hl) in absorption. We limited the redshift range such that the strongest DIBs would fall in the 
optical region. Thus, we selected the highest iV(Hi) DLA galaxies in the redshift range 0.09 < 
z < 0.52. DLAs, by definition, have a large neutral gas column density ( N(Hl) > 2 x 10 20 cm " 2 ). 



However, they are observed to have low reddening (E(B — V) < 0.04) ([Ellison et al.ll2005l ) and 
low metallicity (typically Z ~ 0.1Z & ). The low metallicity and reddening of DLAs suggest that 
the gaseous environments selected by DLA absorption differ from the Galactic ISM in which DIBs 
are observed. As we discuss in § 15.21 the low reddening of DLAs is not merely a bias arising from 
the necessity of low extinction if the background quasar is to be detected at all, since DLAs in 
radio-selected quasars have similarly low reddenings. 

Below, we describe the basic properties of each of the seven DLAs comprising our sample. 
Table [T] summarizes the main parameters, including abundance information. The columns, from 
left to right, list the QSO, absorption redshift of the DLA, column density of neutral hydrogen, 
the zinc abundance, the iron abundance, and the associated references. The Lya line is used to 
measure the DLA redshift. All impact parameters have been converted to a ACDM cosmology 
(H = 70 km s" 1 Mpc" 1 , n m = 0.3, 0\ = 0.7). 

The 

z abs — 0.524 DLA (t^I) toward AO 0235+164 is unique in our sample for many rea sons. 
It is the only DLA known to have DIB absorp tion[j], the A4428 DIB ( Junkkarinen et al. 2004 ) and 
the A5705 and A5780 DIBs jYork et all 120061 ). The A5797, A6284, and A6613 DIBs have limits 
that are not as sensitive due to sky lines. The associated optical galaxy with the smallest impact 



parameter, 6.7 h 7 n kpc, is a late-type s piral and is assumed to be the absor ber ( Rao et al. 



2003 



Burbidge et al.lll996l : lYanny et al.lll989i ). However, IChen Lanzettal (120031 ) argue that many of 
the galaxies in the optical field have the same redshift and may collectively be responsible for the 



lr The DIB-bearing Can absorber reported by Ellison et al. ( 2008h is likely to be a DLA. However, this cannot be 
confirmed without a spectrum of the Lya absorption. 
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metallicities (jTurnshek et al.l 120031 ; lJunkkarinen et al.l 12004!) assume solar abunda nce ratios and 



solar metallicity for the foreground Galactic abs orption. lJunkkarinen et al.l (|2004l ) state that the 



difference between their metallicities and those of Turnshek et al 



DLA absorption. The adopted neutral hydrogen column de nsity is from lJ unkkar inen et al.l (|2004j) 
and is consist ent with iV(Hi) = (4.5 ±0.4) x 10 21 cm" 2 from Turnshek et al.l (|20Q3l ) . The measured 



2003) are likely d ue to different 



realizations of the noise in the X-ray data and/or the variability (jRieke et al.lll976l ) of this QSO. 
Despite the differences in these measurements, this DLA has the highest metallicity in our sample. 
The DLA also has significant reddening with E(B — V) = 0.23, R v = 2.5, and a strong graphitic 
dust feature at 2175 A. 

The z = 0.091 DLA (#2) toward Q0738+313 (01363) is probably a low surface brightness 
galaxy with an impact parameter of < 3.5 hj} kpc ( Turnshek et al.ll2001). This DLA is one of t wo 



found along the QSO sightline of Q0738+313 and was first reported by 



Bao fc Turnshek! (jl99Si ). 



The other DLA, z abs = 0.221 (#3), t oward the QSO Q073 8+313 is probably a dwarf spiral 



with an impact parameter of 20.5 h 7 } kpc ( Turnshek et al.l 12001 



The z abs = 0.518 DLA (#4) toward B2 0827+243 is likely a disturbed spiral galaxy with ex- 
tended gas that p r oduces the observed hydrogen absorption at an impact parameter of 38.2 hy kpc 
(|Rao et al.ll2003l ). iKhare et al.l (J2004J) note that it appears the DLA requires a significant radiation 
field, similar to the radiation within the dense clouds in our Galaxy, to create a low iV(Fei) /iV(Feii) 
upper limit of 10~ 3 . 

The absorbing galaxy giving rise to the z a bs = 0.239 DLA (#5) toward PKS 095 2+179 has not 



been confirmed. However, two candidates lie within < 4.6 kpc (|Rao et al.ll2003l ). The authors 



note that these two galaxies appear to be nearly edge-on and are classified as dwarf low surface 
brightness ga laxies. The DLA appe ars to have multiple Lya components around the central line 
at z = 0.239 (jRao TurnshekH2000l ). The adopted hydrogen column density is measured from the 
central Lya feature. However, the velocity structure suggests a possible clustering. 

The 

z abs — 0.313 DLA (^6) toward PKS 1127—145 is possibly the remnant of a dwa rf low sur 



face b rightness galaxy tidally disturbed by more massive spiral galaxies in the sa me field (jRao et al 



20031 ). The absorber appears to overlap the QSO point sprea d function, thus , 



find an upper limit for the impact parameter of 6.9 h 7 Q kpc. iTurnshek et al.1 (|2003l ) discuss the 



Rao 



et al 



(200 



3) 



difficulty of determining a metallicity for this system due to uncertainties in their X-ray spectrum. 
We have estimated a lower limit of the iron abundance from a VLT/UVES spectrum kindly donated 
for this work by Dr. M. Murphy. The system has a complicated velocity structure. However, fixing 
the Doppler parameter, redshift, and number of Voigt profiles yields a lower limit on the column 
density of log iV(Feli) > 1 5.7 atoms cm -2 . Repeating this analysis using the apparent optical depth 
( Savage &: Sembachlll99ll ) gives a slightly more conservative li mit of log iV(Fe ii) > 15.2 atoms cm -2 , 
which we adopt for the column density of iron. Assuming the lLoddersI (|2003l ) solar abundances, we 
deduce a lower limit of [Fe/H] > — 2. 



The absorbing galaxy responsible for the z abs = 0.395 DLA (#7) along the Q1229-029 sight- 



line is a low surface bri ghtness galaxy with an impact parameter of 8.2 h 7 Q kpc ( Le Brim et al 



19971 : ISteidel et alJll994 ). 



3. Observations and Data Reduction 

Observations of the seven DLAs were obtained with seven facilities toward six QSO sightlines 
between July 2002 and October 2005. The S/N of the QSO spectra range from 5-150 pixel -1 . All 
quoted S /N measurements are calculated in the regions of the expected locations of the redshifted 
A4428, A5780, A5797, A6284, and A6613 DIBs. The journal of observations is presented in Table H 
Cataloged are the QSO, facility and instrument used in the observation, the grating/grism, the 
slit width, the UT date of the observation, the total exposure time in seconds, and the wavelength 
coverage of each spectrum in angstroms. 



3.1. Observations 

QSO spectra covering the z = 0.524, z = 0.239, and z = 0.395 DLAs along the sightlines 
AO 0235+164, PKS 0952+179, and Q1229-020 were obtained with the FORS2 spectrograph on 
the Very Large Telescope (VLT). All resolutions are obtained by measuring unresolved sky emission 
lines. Observations of AO 0235+164 and Q1229-020 use the same grating and the resolutions agree 
within reasonable uncertainties. 

A QSO spectrum covering the z = 0.091 and z = 0.221 DLAs along the Q0738+313 sightline 
was acquired with the DIS spectrograph on the Apache Point Observatory (APO) 3.5 m telescope. 
DIS is configured with a dichroic that splits the light to a blue chip and a red chip at ~5500 A. 
The resolution is measured directly from unresolved sky emission lines. 

A Keck/HIRES spectrum of B2 0827+243 covering the DLA at z = 0.518 was kindly provided 
by Dr. W. Sargent. The resolution is measured by unresolved atmospheric absorption lines (we did 
not have access to sky emission line data for this object). 

QSO spectra covering the z = 0.239 and z = 0.395 DLAs along the sightlines toward PKS 0952+179 
and Q1229-020 were obtained using the ISIS spectrograph on the 4 m William Herschel Telescope 
(WHT). The resolutions were measured directly from unresolved sky emission lines; the uncertain- 
ties are relatively large because of the low S/N in these data. 

A UVES/VLT spectrum covering the z = 0.313 DLA toward PKS 1 127-145 was kindly pro- 



vided by Dr. M. Murphy. The resolution is taken from lDekker et al.1 (120001 ) . We did not have access 
to sky data for this spectrum nor were there any unresolved atmospheric absorption lines, so we 
could not estimate the resolution directly. 



A QSO spectrum covering the z = 0.313 DLA along the sightline toward PKS 1127-145 was 
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obtained with the GMOS spectrograph on the 8.1 m Gemini South telescope. The resolution is 
taken from the online Gemini Science Operations GMOS Instrument Manual. We had limited sky 
data, so we did not estimate the resolution directly. However, due to large equivalent width limits 
of the GMOS data, the effect of uncertainties in resolution should be minimal and not affect results 
in this paper (see Table |3|). 



3.2. Data Reduction 



With the exception of the VLT/UVES and Gemini/GMOS spectra of PKS 1127-145, the data 
were reduced using standard IRAFg packages. The IRAF reduction process involved bias subtrac- 
tion, flat-fielding, spectrum extraction, and wavelength calibration using standard lamps. Once 
wavelength calibrated, each spectrum was continuum fit manually using sfit to achieve the lowest 
residuals in the regions of the DIBs where no detections are observed. Near telluric lines or prob- 
lematic sky subtraction, the continuum was fit using regions bracketing these features. The flux, 
uncertainty, sky (when acquired), and continuum spectra are normalized and then the individual 
spectra are optimally combined (using an algorithm of our own design). For the Gemini/GMOS- 
S spectrum, data reduction was performed using the IRAF Gemini tools in the qmos pa ckage. 
The IRAF task, telluric, was also used on the Gemini/GMOS-S spectrum (jYork et al.l 120061 ) . The 
UVES spectrum was reduced using the standard ESQ pipeline a nd a custom code called the UVES 
Post-Pipeline Echelle Reduction (uves popler. iMurphyibOQO ). 



4. Data Analysis and Spectra 



4.1. Analysis of DIBs 



A modified method of lSchneider et al.l ()1993l ) was developed and used to search for the presence 
of A4428, A5780, A5797, A6284, and A6613 DIB absorption. The Schneider et al. technique is 
optimized for objectively locating unres olved features in s pectr a. However, some DIBs are resolved 
in our spectra. Thus, we modified the ISchneider et al.l (|1993l ) method to be optimized for both 
unresolved and resolved features by combining the DIB's intrinsic full width at half maximum 
(FWHM) with the instru mental spread function (ISF). The procedure is outlined in Appendix A.l. 
The calculation used by ISchneider et al.l (|1993l ) (given in Eq. IA4|) transforms to Eq. IA111 



In addition, we replaced the normalized flux error of lSchneider et al.l (|1993l ) by the residuals 
of the data in pixels where the normalized flux deviates significantly from the continuum (see 
Eqs. IA8l IA91 and lAlOl) . This results in a more conservative detection threshold (equivalent width 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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limit) in the cases of problematic sky subtraction or continuum fits. As an example of this method, 
Fig. [TJa displays the relative flux of Q0738+313 in the region of the expected positions of the A5780 
and A5797 DIBs for the z = 0.091 DLA. Fig. []J contains the sky data. Fig. [lb contains the 
associated 3 a rest-frame detection thresholds for a Gaussian profile with the expected FWHM of 
the redshifted A5780 DIB. Both the A5780 and A5797 DIBs are unresolved. However, if they had 
been resolved they would have yielded different equivalent width limits because the limits depends 
upon the redshifted FWHM for resolved DIBs. The problematic sky subtraction at 6300 A and the 
problematic continuum fit at 6314 A result in conservative limits due to large residuals (see Eqs. lA8| 
IA101 and IA11]) . The region around the A5797 DIB is "clean" in that no residuals are used in the 
equivalent width li mit calculat i ons. The equivalent width limits of the DIBs in the AO 0235+164 



DLA measured by I York et al.l (|2006l ) are less stringent but not inconsistent with those measured 
by this method. 

Possible detections were visually inspected to determine if they were a DIB or perhaps an- 
other feature or sky line residual. Checking each candidate feature is essential because of a small 
uncertainty in the rest wavelength of each DIB as well as a small uncertainty in the redshift of DIB 
absorbing gas relative to the Lya determined redshift of the DLA. For most DIBs, we were able to 
measure only equivalent width limits. 

To quantify confidence levels in the measured detection thresholds determined with our mod- 
ified method, we have estimated the 1 a uncertainties in the measured equivalent width limits. 
There is some uncertainty in the resolution (ISF) of each spectrum. We estimate the uncertainty in 
the resolution as the standard deviation of the A/FWHM ratios of unresolved sky absorption lines or 
emission lines, where the A are the line centers. If sky data are not available, we excluded the uncer- 
tainty in resolution. In addition, there is uncerta i nty in the accuracy of the continuum fits, which we 



estimated using the technique of lSembach et al.l (|199ll ). Of importanc e, is the fact that DIBs hav e 



measured uncertainties in their rest-frame wavelengths and FWHMs (IJenniskens Sz Desertl 11994). 
The uncertainty in the equivalent width limit due to uncertainty in the wavelength is the standard 
deviation of the individual equivalent width limits computed over the range of the redshifted DIB 
wavelength uncertainty. The full explanation of "error" propagation to obtain uncertainties in the 
equivalent width limits is given in Appendix A. 2 (see Eq. IA12j) . 

Contained within Fig. [1] is an example of our equivalent width limit analysis including infor- 
mation on the effects of uncertainties in the wavelength of the band center. If a redshifted DIB 
wavelength is near a problematic region then an uncertainty in rest wavelength can introduce a large 
uncertainty in a measured equivalent width limit. This is a more substantial issue for the A4428 
DIB, which has the largest uncertainty in its rest-frame wavelength (~ 1.4 A). For the majority of 
the equivalent width limits, the largest uncertainty is due to the continuum fit. 



Only three DIBs have been detected in DLAs, the A4428 DIB reported in I Junkkarinen et al 

200- 

( 



(2004), and the A5705 and A5780 DIBs reported in a companion paper of this work by 



York et al 



20061 ) . All three arise in the z = 0.524 DLA toward AO 0235+164. The equivalent widths were 
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calculated by sum ming the individual equivalent widths of each pixel. Both of the detections from 



York et al.l (|2006l ) were at least one resolution element away from sky lines. The undetected A5797 
and A6284 DIBs toward the z = 0.524 AO 0235+164 DLA were located within strong sky lines. 



York et al.l (J2006J) estimated the equivalent width limits for these DIBs using sy nthetic abs o rption 



features with varying FWHMs. The undetected A6613 DIB is near a sky line; lYork et al.l ((2006) 
estimated its equivalent width limit directly from the signal-to-noise, dispersion of the chip, and 
number of pixels. 



4.2. Measurements and Spectra 



The measured equivalent widths, equivalent width limits (3 a), and their 1 a uncertainties are 
listed in Table [3l Tabulated are the DLA (by number), the corresponding QSO, the absorption 
redshift of the DLA, the facility and instrument, and the rest-frame equivalent widths or limits, 
with uncertainties (mA) of each of the DIBs if they were observable. Spectra from which the most 
stringent limits were obtained are shown in Figs. [2HcfR Each panel (a)-(e) displays the region 
around a redshifted DIB. The upper sub-panels display the normalized flux (histogram) with the 
expected positions of the DIBs (marked by ticks) based upon the Lya redshift. With the exceptions 
of Figs. [5] and d the center sub-panels display the uncertainty spectra of the associated normalized 
fluxes, and the lower sub-panels display the sky counts normalized by the continuum. Sky data 
were unavailable for B2 0827+243 at z = 0.518 and PKS 1127-145 at z = 0.313. 

The smooth thin curves through the data are the expected observed DIB absorption profiles 
based on the measured TV (Hi) for the PL As, where the band centers and FWHM of each DIB 
are taken from iJenniskens Desert] (| 19941 ) . These models are not computed for the A4428 and 
A6613 DIBs since these DIBs have no published iV(Hi) relationships. In the following sections we 
will discuss limits on reddening and gas-to-dust ratios. Thus, we also illustrate (thick curves) the 
expected observed DIB absorption profiles assu ming an E(B — V) = 0.04, the upper limit for high 
redshift DLAs assuming SMC-like extinction dEllison et al.l hood) . F or the AO 0235+164 DLA, 
we adopted the measured E(B — V) = 0.23 (jJunkkarinen et al.ll2004l )). The computations of the 
reddening models are discussed in § 15.21 



The equivalent widths of the DIBs esti mated from the Galactic DIB-iV(Hl) and Galactic 
BIB-E(B - V) models (see lWelty et alJbood ) are presented in Table H Columns 1-2 list the DLA 
number and the QSO with associated DLA redshift. Columns 4-8 provide the observed equivalent 
widths and equivalent width limits for the A4428, A5780, A5797, A6284, and A6613 DIBs. Also listed 
are the predicted equivalent widths of the DIBs, where EW[iV(HI)] denotes the Galactic iV(Hi) 
scaling and EW[E(B — V)] denotes the reddening scaling in mA. E(B — V)n m is the reddening 



3 The VLT/UVES spectrum is not shown because only the A4428 DIB is covered. The UVES spectrum has a high 
resolution that makes identifying the broad A4428 DIB very difficult. 
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upper limit for each DLA based upon the observed DIB equivalent width limits and the reddening 
relation (see § 15 . 2[> . 



5. Results and Discussion 

In this section, we examine the DIB strengths in DLAs compared to the Galactic DIB-A^(Hi) 
relation and the Galactic E(B — V) relation. The former provides information on the gas content 
and the latter provides information on the reddening (an indirect indicator of dust content). Our 
observations allow us to estimate lower limits on the gas-to-dust ratios, N{R\) / E{B — V), of the 
DLAs in our sample. 



5.1. Gas Content 



The widely observed Galactic DIB-A r (Hi) relation describes the correlation of the equivalent 
width of the A5780 , A5797, and A6284 DIBs with th e column dens i ty of neutral hydrogen along the 
same line of sight (jHerbigi Il995l : IWelty et al.l 12006 ) . IWelty et al.1 (120061 ) extend the Galactic work 
by including Magellanic Cloud sightlines. 

Plotted in Fig. M are the DIB-N(Hi) relations for the A5780 (panel a), A5797 (panel b), 
and A6284 DIBs (panel c), where logiV(Hi) [cm -2 ] is plotted against the logarithm of the DIB 
equivalent widths [mA]. Also plotted are the (1 a) error bars for measured values or downward 
arrows representing upper limits. The vertical error bars for the A5780 DIB detection in the 
AO 0235+164 DLA are s maller than the poin t size. The solid lines represent the weighted best-fits 
to the Galactic data from lwelty et all ( 2006 ). The equivalent widths predicted at a given iV(Hi), 
the EW[iV(Hi)] listed in TableHJ are computed using this best-fit line and the A r (Hi) of each DLA 
(Figs. [2HE] contain these models superimposed on the data). The dotted lines in Fig. [9] roughly 
enclose the regions containing the Galactic data. The dashed lines roughly enclose the regions 
containing the LMC data, and the dot-dashed lines roughly enclose the regions containing the 
SMC data. 

The DIBs in the DLAs whose equivalent width limits in Fig. [9] lie below the Galactic best-fit 
lines are robust enough that we should have detected them if they followed the same dependence 
on iV(Hi) as Galactic sightlines. In several DLAs, the limits are inconsistent with the expectations 
from the Galactic sightlines; the DIB strengths are weaker than expected. The DIB limits are 
consistent with the strengths of DIBs in the LMC or SMC. However, higher signal-to-noise data 
are required to determine whether the DIB strengths are actually consistent with or are even weaker 
than those observed in the LMC and SMC. The A6284 DIB, panel (c) of Fig. [9j provides the most 
stringent evidence that the DIB strengths in four DLAs are not consistent with DIB strengths 
in the Milky Way. The four DLAs toward AO 0235+164 (z = 0.524), Q0738+313 (z = 0.091), 
PKS 0952+179 (z = 0.239), and PKS 1127-145 (z = 0.313) are underabundant in their A6284 DIB 
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strengths by factors of 4-10 times expected from the Galactic DIB-iV(Hi) relation. 



A number of environmental factors such as reddening, metallicity, and ionizing radiation may 
be responsible for DLAs not following the Galactic DIB-iV(Hi) relation. Whether it is one or a 
combination of these factors, the data suggest that the environments probed by DLA sightlines 
differ from those within the Milky Way in which DIBs are present. The environments that give rise 
to the DIBs are likely very localized. Galactic properties vary even on small scales, which is why 



Cox et al.l (|2007l ) find varying differences and DIB strengths along different sightlines within the 



larger confines of the SMC. QSO sightlines probe relatively small transverse spatial scales in DLA 
galaxies. Therefore, our results do not eliminate the possibility that DIBs follow the Galactic DIB 
iV(Hi) relationships elsewhere in these galaxies. We are less inclined to suggest that our results 
indicate redshift evolution of organics, since infrared emission from PAHs has been o bserved in the 
extremely dusty environments of ULIRGs as high as redshift z ~ 2 (|Yan et al.ll2005l ). 



5.2. Reddening 



There has bee n a long history of investigating DIB dependence on color excess ([Merrill fc Wilson 



1938 



Herbiglll993l ). Whereas Galactic sight-lines along the disk (low Gala ctic latitudes) ha v e typ 



ical reddening values of 0.1 to 1.0, DLAs typically have lower reddening. iMurphv 8z Liskd (|2004! ) 
use Sloan Digital Sky Survey (SDSS) QSO spectra to esti mate E(B — V) < .02 mag (3 a) along 
DLA sightlines at z ~ 3, assuming SMC-like extinction. lEllison et al.l (|2005l ) use radio selected 
QSOs from the Complete Optical and Radio Absorption Line System (CORALS) survey to esti- 
mate reddening along DLA sightlines. The CORALS survey avoids the potential problem of bias 
from optical luminosity selected samples that may inhibit the detection of more reddened QSOs. 
They find E(B - V) < 0.04 (3 a) assuming SMC-like extinct ion for 1.9 < z < 3.5 DLAs. Using 
a sample of 0.8 < z < 1.3 DLAs selected by Call absorption, IWild &: Hewettl (120051 ) estimate an 
average reddening of E(B — V) = 0.06 assuming LMC and SMC extinction curves. 

The correlations between reddening and the A5780, A5797, and A6284 DIB str engths are tight 
among Galactic and Magellanic Cloud sightlines, with best fits (jWelty et al.l l2006). 



logEW[E(B - U)]a5780 = 0.99 log E(B — V) + 2.65, 



(1) 



and 



logEW[E(B - IOU5797 = 0.99 log E(B - V) + 2.26, 



logEW[£(£ - y)] A6 284 = 0.80 log E(B -V) + 3.08. 



(2) 



(3) 



The rms scatter of the measured equivalent widths about these relationships are less than 0.15 dex; 
in other words, the A5780, A5797, and A6284 DIBs fairly e qually obey th e BIB -E(B - V) relation 
in the Galaxy. Combining Galactic and extragalactic data, lEllison et al.1 (J2008J) find a A5780 DIB- 
E{B — V) relation with a slightly larger slope, 1.27, and intercept, 2.77. Compared to lWelty et al 
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(|2QQg), the Ellison et al.l (j2008h results yield a 10% difference in logEW[E(B - F)] A 5780 for unit 
reddening. 

Correlations with the A4428 DIB equivalent width are notoriously problematic to compute 
because of the difficulty in measuring the DIB's broad width. Previous wo r k has noted a correla- 
tion of the A6613 D IB equivalent width with E(B - V) (IMegier et allhood : iThorburn et al.ll2003l : 
Weselak et ahlboOlh . In the case of both the A4428 and A6613 DIBs, we employ the reddening 
relations provided by T.P. Snow (unpublished, private communication). Using the best-fit for 75 
Galactic sightlines for A4428 and 123 Galactic sightlines for A6613, 



EW[E(B - V)]xu28 = 2093.99^(5 - V), 



(4) 



and 



EW[E(B - V)] 



A6613 



2nmE{B - v). 



(5) 



The scatter about these relations is relatively large when compared to the other DIBs studied in 
this work. The A4428 DIE>-E(B — V) relation has a Pearson correlation coefficient of 0.825, while 
the A6613 DIB-E(B — V) relation has a Pearson correlation coefficient of 0.831 (unpublished, T.P. 
Snow, private communication). 



Assuming the upper limit measured for high redshift DLAs of E(B—V) = 0.04 from lEllison et alJ 



(|2005l ). we computed the expected DIB equivalent widths using Eqs. [H El El HI and El The 
only exception is the AO 0235+164 DLA, which has a measured reddening of E{B — V) = 0.23 
(IJunkkarinen et al.1 l2004h . The expected DIB equivalent widths are listed in Table HI Assuming 
these reddening values, the expected DIB equivalent widths are systematically smaller than those 
predicted by the DIB-iV(Hi) relationship. In Figs.EHSwe show the expected DIB profiles (thick 
curves). With the exception of the A6284 DIB limit measured from the spectrum of the DLA 
toward AO 0235+164 and the A4428 DIB limit measured from the spectrum of the DLA toward 
PKS 1127-145, the expected DIB equivalent widths for adopted reddening values are smaller than 
our measured equivalent width limits; thus, if these reddening values are appropriate for interme- 
diate redshift DLAs, we would not have detected these DIBs. As discussed in § [U the A6284 DIB 
limit for the AO 0235+164 system is problematic due to a large atmospheric feature and the A4428 
DIB limit for the PKS 1127-145 system is unreliable due to difficulties in continuum fitting the 
broad DIB within the high resolution VLT/UVES data. 

Assuming the Galactic T)IB~E(B — V) relations are valid for DIBs in DLAs, we used our 
measured equivalent width limits for the A4428, A5780, A5797, A6284, and A6613 DIBs to compute 
upper limits on the reddening for the DLAs in our sample d irectly from t he we ighted slopes of 
the BIB-E(B - V) relationships (Eqs. CO El El H and [5]) from IWeltv et al.1 (|20od ) and T.P. Snow 
(unpublished, private communication). The reddening limits, E(B — V)ii m , are listed in Tabled] 
for each DIB. 



In Fig. [TUl we illustrate the sensitivity of this method for E(B — V) values of 0.02, 0.04, 0.10, 
0.20, 0.40, and 1.0, respectively. Plotted are log(EW Limit /EW[E{B - V)}) for each DLA for the 
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A4428, A5780, A5797, A6284, and A6613 DIBs (panels a, b, c, d, e, respectively). For the condition 
log(EW Limit /EW[£?(5 — V)]) < 0, the reddening predicted equivalent width for a given DIB is 
greater than or equal to the upper limits afforded by our data. A value of log(EW Limit /E W [E (B — 
V)]) = corresponds to our computed reddening upper limit, E(B — F)ii m , at which the DIB 
should just become detectable in our data. The vertical error bars define the 1 a uncertainties in 
log(EW Limit/EW[£?(5 — V)]) evaluated at our computed E{B — V)\\ m . The uncertainties are cal- 
culated using standard error propagation taking i nto account the unc ertainties in both the W Limit 
and the slope of the ~EW[E(B — V)] relations from lWelty et al.l (|2006l ) and T.P. Snow (unpublished, 
private communication). The A4428 and A578 DIBs for the PL A tow a rd AO 0235+164 are left 



blank because they are confirmed det ections ([J unkkari nen et al. 
known reddening of E{B — V) = 0.23 (j Junkkarinen et al.l 12004 ) . 



2004 



York et al.l l2006h with a 



Note that in Fig. [T0| the A6613 DIB-E(B — V) relation is the least constraining of the DIBs in 
determining the upper reddening limit, E(B — V)n m . Furthermore, due to the scatter in the A4428 
and A6613 DIB-E(B — V) relations, the errors in the best-fit slope in Eqs. [Hand [5] are relatively 
large which results in the large errors displayed in Fig. [TUb andllOb. For these reasons, we do not 
adopt upper reddening limits, E{B — V)^, from the A4428 or A6613 DIBs. 

Our final adopted reddening limit for a given DLA is determined from either the A5780 or 
A6284 DIBs. The adopted limits are listed in Table [5j Two of the DLAs have limits of E(B — V) < 
0.05 (Q0738+313 z = 0.091 and PKS 0952+179), both from the A6284 DIB. The DLA toward 
B2 0827+243 has an upper limit to the reddening of E(B - V) < 0.07 from the A5780 DIB (the 
A6284 DIB is not covered). The A6284 DIB provides upper reddening limits for the remaining 
DLAs toward Q1229-020 {E(B — V) < 0.08), Q0738+313 z = 0.221 (E(B — V) < 0.14), and 
PKS 1127-145 (E(B - V) < 0.21). If we were to apply constraints from the A4428 DIB, we would 
infer upper reddening limits of 0.07, 0.04, and 0.03 for the z = 0.221 DLA toward Q0738+313, and 
the DLAs toward B2 0827+243 and PKS 1127-145, respectively. 



For the DLA at z = 0.524 toward AO 0235+164, there are two detections from lYork et al 



(|2006j) (the A5705 and A5780 DIBs) and one detection from l.lunkkarinen et al.1 (|2004]) (the A4428 
DIB). The equivalent widths measure d from the A4428 and A5 780 DIBs are larger than expected 
from the measured E(B - V) = 0.23 (j Junkkarinen et al.ll20o3 ). The A5797 and A6613 DIBs give 
limits consistent with this, but the A6284 DIB suggests an E{B — V) < 0.06 by this technique. 
Essentially, we should h ave discovered the A6284 DIB assuming our technique is correct because 
Junkkarinen et al.1 (J200J) measures an E(B — V) = 0.23 ± 0.01. 



A plausible reason for our non-detection of the A6284 DIB in the AO 0235+164 DLA is that the 
expected position of the A62 84 DIB at z = 0.52 4 resides directly in a broad atmospheric absorption 
band (see panel b in Fig. [2|) . I York et al.l (|2006l ) attempted to overcome this by modeling the A6284 
DIB with various equivalent widths and convolving it with the atmospheric band. The limit they 
achieve (see Table [3|) is based on the minimum equivalent width needed to separate the DIB out of 
the atmospheric band. This is a robust way to measure the limit. However, there is some inherent 
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difficulty in this method because it has to be assumed a priori that the atmospheric band is not 
contaminated by DIB absorption. A second potential reason that we do not detect the A6284 DIB 
in the AO 0235+164 DLA, given its high reddening, is that this DIB may not follow the Galactic 
DIB-E(B — V) relation in DLAs. Indeed, in the Magellanic Clouds, A6284 DIBs are weaker than 
their Galactic coun terparts by a factor of two relative to the Galactic A6284 DIB— E(B — V) relation 
(jWelty et al.ll2006l ). Given that the A6284 DIBs in four of our DLAs are constrained to be 4-10 
times weaker than the Galactic A6284 DIB-iV(Hi) relation, it is reasonable that DIB relations in 
DLAs could be more Magellanic Cloud-like than Galactic-like. 

Our observations have provided reddening constraints of 0.05 < E(B — V) < .08 fo r four 



DLAs at 0.09 < z < 0.52. This extends to lower redshift t he work of lEllison et al.1 (|2005l ). who 
report E(B — V) < 0.04 for 1.9 < z < 3.5 DLAs, and of iMurphv fc Liskel ((20041 ). who report 
E{B — V) < 0.02 for z ~ 3 DLAs (all limits are 3 a). Reddening is an indirect measure of dust 
content. The slope of any T)IB-E(B — V) relation li kely has a strong depende nce on the nature of 
the dust (particle size, abundance, and composition). lJunkkarinen et al.1 (|2004l ) argue that the dust 
in the AO 0235+164 DLA is more Galactic-like than Magellanic-like, but with fewer small particles. 
Because the organics responsible for the DIBs may depend on dust for formation and/or survival, 
the abundance and nature of the dust in the DLA toward AO 0235+164 might be responsible for 
the presence of DIBs observed in this DLA. 



5.3. Gas-to-Dust Ratios 



Assuming that our method of estimating E(B — V)ii m realistically reflects the upper limit on 
reddening in our DLAs, we can estimate lower limits in their gas-to-dust ratios, N(Hi)/E(B — V). 
We present the gas-to-dust ratio lower limits in Table G2 they are computed using the measured 
iV(Hi) from the literature and our adopted upper E{B — V) limits. The lower limits range from 
2.9 x 10 21 to 42 x 10 21 cm -2 mag" 1 . Note that the gas-to- dust ratio 19.2 x 10 21 cm " 2 mag -1 for 
the z = 0.524 DLA toward AO 0235+164 was measured by I Junkkarinen et al.l |200J). 



Bouchet et al.1 ( 19851 ) measure lower a nd upper value s, 37 x 10 21 and 52 x 10 21 cm 2 mag 



-1 



respectively, for SMC gas-to-dust ratios. iGordon et al.l (120031 ) determined a gas-to-dust ratio of 
19.2 x 10 21 cm" 2 mag" 1 for the LMC (from their LMC-2 data). T he average LMC s ample of 

A linear fit by Icox et al.l (j2006bl ) to their 



Gordon et al.1 ((20031 ) yields 11.1 x 10 21 cm" 2 mag" 1 



LMC data g i yes a similar gas-to-dust ratio of 14.3 x 10 cm mag . F or the Milky Way , 
Bohlin et al.l (|1978l ) find a gas-to-dust ratio of 4.8 x 10 21 cm -2 mag -1 ; whereas. I Cox et al.l ((2006b) 



find a ratio of 4.03 x 10 21 cm 2 mag 1 from the fit to their Galactic data. 

Although the SMC, LMC, and Milky Way appear to have distinct ranges of gas-to-dust ratios, 
in reality, there is a continuum that likely reflects the variation of local environments in each. 
Gordon et al.1 (120031 ) measure E(B — V) and A r (Hi) toward four stars in the bar of the SMC and 
one star in the wing of the SMC. From these data, the gas-to-dust ratios in the SMC bar span 
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th e range 17-51 x 10 21 cm 2 mag 



T hese va lues overlap with both the LMC gas -to-dust ratios 



of iGordon et al.l (|2003l ) and ICox et al.l (|2006bl ) and the SMC gas-to-dust ratios of IBouchet et al 



1985JV The SMC wing data yield 15.2 x 10 



21 cm 2 



mag 1 , which falls very near the ICox et al 



(j2006bh LMC fit. 



Cox et al.l (120071 ) searched for DIBs toward six lines of sight in the bar and wing of the SMC. 



The sightline AzV456, which is in the wing of the SMC, is the only sightline in which they detected 
DIBs. From their E(B — V) and 7V(Hi) data (Table 9 on-line material), the gas-to-dust ratio in 
the SMC wing is ~ 7 times lower (~ 7.4 x 10 21 cm -2 mag -1 ) than the gas- to-dust ratio in the 
SMC bar (~ 52 .4 x 10 21 cm" 2 mag -1 ). Th is SMC wing value li es be tween thebohlin et al.l (|l978h 



and lCox et al.1 (|2006bJ) Milky Way and the IGordon et al.l (120031 ) and ICox et al 



(200 6]]) LMC gas - 



to-dust ratios. The SMC bar value lies just above t he upper range o f the IBouchet et al.1 (jl985l ) 
SMC gas-to-dust ratio. It is of interest to note that ICox et al.l (120071 ) detected DIBs only in the 
sightline with the smallest gas-to-dust ratio in their sample and that this ratio approaches the range 
observed for the Milky Way. 



Cox et al.l (|2007l ) argue that the SMC wing is quiescent w hile the star f ormat ion regions of 
the SMC bar are turbulent and exposed to larger UV fluxes. IWeltv et al.l (|2006l ) calculate an 
increase of ~ 28-83 times the average interstellar radiation field (ISRF) for sightlines in the SMC 
bar near star-forming Hn regions; whereas, the SMC wing is more similar to the Galaxy with a 
radiation level of ~ 0.6 times the average ISRF. The variation in the gas-to-dust ratios within the 
SMC are lik ely a reflection o f the balance of dust formation processes, such as accretion, with dust 
destruction ( Cox et al. 2007h. Thus , the dust, and possibly the organics, may be destroyed in the 



SMC bar explaining why I Cox et al.l (|2007l ) could not detect DIBs in their sample of SMC bar stars. 



In Fig. [TTJwe plot log iV(Hi) versus log E(B — V)n m for the DLAs in our sample. The 7V(Hi) 
data points are taken from Tableland the E(B — V)n m are taken from TableEl The log iV(Hi) error 
bars are 1 a. The leftward arrows represent reddening upper limits for our sample. For Q0738+313 
at z = 0.221, B2 0827+243, and PKS 1127-145, the A4428 DIB provides more stringent limits on 
E(B — V)n m . However, they should be considered less robust due to the large errors as shown in 
Fig. HOl panels a and e. 

Although the gas-to-dust ratios in the SMC, LMC, and Milky Way exhibit some overlap, 
in general the Milky Way IS M values tend to be lower than those of the LMC which tend to be 



lower than those of the SMC ( Bohlin et al 



19781 : 



Bouchct et al. 



1985 



Gordon et al 



2003; 



Cox et al 



Bouchet et al 



2006b). In Fig.CCQ the dot-dashed lines illustrate the upper and lower SMC values of 
(|l985l ). The middle three lines in Fig. [TTJ represent gas-to-dust ratios for the LMC, and the bottom 
two lines re present those for the Galax y. The long-dashe d line gives the LMC-2 gas-to-dust ratio 



derived from I Gordon et al.l (|2003l ) . The 



line. The short-dashed line is the 



Gordon et al 



dashed MW line is the iBohlin et al 
ratio. 



Cox et al.l (l2006bh ratio for the LMC is given by the dot-dot 
( 20031 ) ratio from their average LMC sample. The 



dl978|) ratio, and the dotted MW line is the lCox et al.l ((2006a) 
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As shown in Fig. [TT1 the DLA toward PKS 0952+179 has a gas-to-dust ratio that is consistent 
with or greater than the SMC. Two of the DLAs in our sample, Q0738+313 (z = 0.091) and 
PKS 1127-145, have ratios consistent with or greater than the LMC gas-to-dust ratios. Two, 
Q0738+313 (z = 0.221) and Q1229-020, are consistent with or greater than the Galactic ratios. 
The ratio for B2 0827+243, which has the lowest TV (Hi) in our sample, is not well constrained. 
The measured ratio for the z = 0.524 D LA toward AO 0235+164 is iV(Hl) /E(B - V) = 19.2 x 10 21 



cm 2 mag 1 ( Junkkarinen et al .1120041 ). which is consistent with the LMC-2 ra tio of I Gordon et al 



(J2003J) and four times larger than the Galac tic ratio from iBohlin et al.1 (|1978l ). This DLA has a 



very large iV(Hi), and the largest E(B — V) (jJunkkarinen et al 



2004) of the DLAs presented here. 



Three DLAs are constrained to have slightly larger gas-to-dust ratios than the z = 0.524 DLA 
toward AO 0235+164, the z = 0.091 DLA toward Q0738+313 (~ 1.5 times greater), z = 0.239 
DLA toward PKS 0952+179 (~ 2.2 times greater), and the z = 0.313 DLA toward PKS 1127-145 
(~ 1.2 times greater). However, for the DLA toward PKS 1127-145, if we apply E(B — V)y im from 
the A4428 DIB, the gas-to-dust ratio is ~ 7 times greater than that of the AO 0235+164 DLA and 
would be more consistent with gas-to-dust ratios greater than measured in the SMC. Again, we 
caution that this more stringent ratio is not adopted due to large errors from the scatter in the 
A4428 BIB-E(B - V) relation (see Eq.HJ). 

The sightlines toward DLAs are certainly gas rich, which they share with Galactic sightlines 
that contain DIBs. However, their environments, as probed by the background QSO, may have large 
variations in E(B — V), metallicity, and radiation. Understanding the magnitudes and properties of 
the dust, radiation, Hi and H2 content, and metallicity in DLA samples holds promise of revealing 
whether the environments of DLAs are conducive for the organics that give rise to DIBs. Reddening 
is seemingly low in the sightlines probing the galaxies in this sample, which may play a pivotal 
role in inhibiting DIB strengths. The strengths of DIBs in DLAs do not follow the Galactic iV(Hi) 
relation (see Fig. |9|). However, the A5780 detection in the AO 0235+164 DLA is co nsistent with 
expe ctations from the DIB-E(B — V) relation for Galactic and extragalactic points (lEllison et al 



20081 ). Thus, gas-to-dust ratio may not be as good a predictor as reddening for the presence of 



DIBs in DLAs. This would suggest that DIBs might be present in DLAs with a wide range of 
gas-to-dust ratios, but that they lie within the right-hand region of the log iV(Hi)-log E{B — V) 
plane of Fig. [TT1 



Ellison et al.l (|2005l ) conclude that their sample of 1.9 < z < 3.5 high redshift DLAs does not 



follow a Galactic gas-to-dust relation, which is consistent with what we find for four of the DLAs 
in this work (AO 0235+164 at z = 0.524, Q0738+313 at z = 0.091, PKS 09 52+179 at z = 0.239, 
and PKS 1127-145 at z = 0.313). However, as found by lEllison et al.1 (|2005l ). we cannot conclude 
whether the gas-to-dust ratios in our lower redshift DLAs are consistent with the SMC gas-to-dust 
ratios. The remaining three DLAs (Q0738+313 at z = 0.221, B2 0827+243, a nd Q1229-020) do n ot 
have sufficient reddening limits to compare with the high redshift sample of lEllison et al.l (120051 ). 
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5.4. Metallicity versus Reddening 



Cox et al.l ([20071 ) give evidence that metallicity plays an important role in DIB strength since 
a high carbon abundance is necessary for the creation of the organic molecules and the dust grains 
on which they may be formed. In Table [H we list the zinc and iron abundances from the literature 
for our DLA sample. Zinc is a good metallicity indicator because it traces iron-group abundances 
and does not readily deplete on dust. Unfortunately, only upper limits on zinc abu ndances have 
been measured, with the exception of the Q1229-020 DLA, where [Zn/H] = —0.47 (jBoisse et al 
19981 ). Thus, we cannot address any trends of DIB strengths with metallicity. 



The Q1229-020 DLA has a metallicity roughly 0.5 dex above the average metallicity of DLAs, 
and its metallicity is consistent with that of the AO 0235+164 DLA. Since it has a relatively 
stringent upper limit on the reddening, E(B — V)u m < 0.08, much below the E{B — V) = 0.23 
of the AO 0235+164 DLA, it may be a promising candidate for contributing to an understanding 
of the role of metallicity in determining DIB strengths in DLAs. If metallicity plays a role, it is 
expected that DIB strengths in DLAs would to some degree scale with iV(Hl). The JV(Hi) of the 
Q1229-020 DLA is a full dex below that of the AO 0235+164 DLA. Thus, that our upper limit 
on the equivalent width of the A5780 DIB in the Q1229-020 DLA is a factor of 1.7 below that 
detected for the AO 0235+164 DLA is not constraining. A very deep spectrum of the DIBs in the 
Q1229-020 DLA would be very interesting. 

Selecting DLAs by metal absorption would represent another approach to se arching for DIBs in 
external galaxies, since these DLAs may represent a more reddened population. IWild et all ([2006) 
propose that Can absorbers may represent an intermediate link between the quiescent, metal-poor 
and dust-poor DLAs and the intermediate redshift, star- forming, and metal -rich Lyman break 
galaxies with typical reddening of E(B — V) = 0.15 to 0.20. Ellison et al.l (|2008l ) searched for DIBs 
in nine 0.07 < z < 0.55 Can-selected absorbers, and detected the A5780 DIB in only one. At this 
early juncture, the Can selection success rate is comparable to that of Hi-selected DLAs. Mgll 



absor bing DLAs may also represent a slightly more reddened population of DLAs ([Menard et al 
2008h . 



6. Conclusions 



In this paper we employ a generalized method of the ISchneider et al .1 (|l993l ) technique to find 
lines and determine conservative equivalent width limits for DIBs in seven DLAs along with an 
assessment of uncertainties in these limits. We find: 

(1) The A6284 DIB in four of the DLAs in our sample have equivalent width upper limits that 
are 4-10 times lower then expected for similar A^(Hi) relative to Galactic sight lines. These limits 
are not inconsistent with the A6284 DIB strengths found in the LMC and SMC. 



(2) Assuming the A5780 and A6284 BIB-E(B - V) relations hold for DLAs, as it does for 
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Galactic and Magellanic Cloud sightlines, we estimated upper limits on the reddening for our 
sample. In four of our DLAs we estimate E(B — V) < 0.08 with two DLAs having an E(B — V) < 
0.05. These results are consistent with hig h redshift DLA sample s, which have E(B — V) < 0.04 
([Ellison et alJbood ) and E(B - V) < 0.02 (jMurphv fc Liskelbooi l- 



(3) Applying the E(B — V) limits, one of our DLAs is consistent with having the same or larger 
gas-to-dust ratio as the SMC. Two of our DLAs are consistent with having gas-to-dust ratios at 
least as large as sightlines in the LMC. Three of our DLAs have less stringent limits that give lower 
limit gas-to-dust ratios consistent with Galactic, LMC, or SMC sightlines. The AO 0235+164 DLA 
has a measured E(B — V) and iV( Hl) that puts its gas-to-d ust fraction on the high end of the LMC 
sightlines as stated in other work (jjunkkarinen et al.l 12004 ). The E(B — V)n m < 0.06 constrained 



from t he A6284 DIB is inconsistent with the known E(B — V) = 0.23 measured in 



Junkkarinen et al 



(j200J ). We shoul d have been able to detect the A6284 DIB in the AO 0235+164 DLA given the 
limit adopted by York et al. (j20od ). Three possibilities are that our method of determining the 
E{B — V)ii m does not apply to DLAs, our limits do not adequately take into account the large 
atmospheric absorption band (see Fig. [2f)), or that conditions are not favorable to the formation or 
survival of this DIB carrier. 

It is interesting to speculate whe ther ionizati o n con ditions may be an important factor in 
inhibiting or enhancing DIB strengths. IWelty et al.1 (120061 ) test the ionization effects of DIBs along 
Galactic, LM C, and SMC line s of sight. However, they do not find any significant trends. On 
the contrary, ICox et al.1 (|2007l ) measure the UV radiation field along lines of sight toward the 
SMC and find UV radiation is an important environmental factor in DIB strengths. Laboratory 
spectroscopists claim that the DIBs may be due to parti ally ionized PAHs because they produce 
a wealth of absorption features in the optical spectrum (jSnowl l200ll ). If this is the case then a 



significant UV radiation may be required. However, UV radiation that is too high will dissociate 
the molecules. In this work, we are unable to explore the affects of ionization conditions; knowledge 
of ionization conditions in specific DLAs is difficult to obtain. 

It is also interesting to speculate whether metallicity may be an important factor in inhibiting 
or enhancing DIB strengths. However, since we have robust metallicity measurements for only two 
of the DLAs in our sample, we cannot directly address the affects of metallicity. We do point out 
that obtaining deep spectra of DIBs in DLAs with known metallicity could be a fruitful future 
research direction, especially if reddening were also known. A direct comparison between the low 
metallicity Q1229-020 DLA and the high metallicity AO 0235+164 DLA might be fruitful as a first 
examination of the affects of metallicity in determining what governs DIB strengths in DLAs. 



6.1. Implications of our Results 



Our results imply that reddening is a more crucial indicator of DIB strengths than is Hi content 
for DLAs; the low reddening in DLA selected galaxies inhibits the presence of DIBs. However, our 
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sample is small and additional observations are required to ascertain the strength of this statement 
(metallicity and ionization conditions likely play a role as well). 

The weakness of DIB strengths in our sample hints that the environments of high iV(Hl) DLA- 
selected galaxies may be less suitable to create and/or sustain the organic molecules than those of 
the Galaxy. Not only is the immediate solar environment beneficial for sustaining life, but it may 
be that the Galaxy is a more hospitable location for the survival of organic molecules that may 
have been the precursors to biology on Earth. If these molecules are important as precursors to life 
in the universe, charting their presence to high redshift places constraints on how long ago and in 
which environments life could have potentially formed. The presence of DIBs in the AO 0235+164 
DLA demonstrates that organic molecu les existed in at least one DLA selected environment at a 



redshift of z ~ 0.5, or some 5 Gyrs ago (jYork et al.l l2006). 



Due to their weakness relative to DIBs in the Galaxy, observing DIBs in the general population 
of high redshift galaxies (as opposed to ULIRGs and star bursting galaxies) remains a challenge. 
Selecting galaxies by DLA absorption may yet prove to be a lucrative method for detecting DIBs 
at high redshifts, provided the detection sensitivity can be increased. However, most known DLAs 
reside at high redshift where the DIBs move into the near-IR, where high sensitivity spectroscopy 
is time intensive. 
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A. Appendix 

A.l. Measuring Equivalent Width Detection Thresholds 



A straight forward method for determining equivalent width limits was presented bv lLanzetta et al 



(jl987l ). The uncertainty (or error) spectrum is used to compute the equivalent width limit in each 
individual pixel, i, 

EM (Al) 



0",, 



D, 



I c {\ 



where Di is the pixel dispersion in angstroms, E(Xi) is the uncertainty in the flux at Aj, and I c (Aj) is 
the estimated continuum flux at Aj. The 1 a equivalent width detection threshold (limit) centered 
at Aj is then obtained by summing the individual equivalent width limits of adjacent pixels over a 
selected aperture, 



.12 



3=31 



1/2 



(A2) 



where j\ = i — m/2, ji = i + m/2, and where m, the aperture size, is an even number of pixels 
over which the absorption feature is anticipated. For unresolved features, m can be taken to be 
roughly two resolution elements (i.e., roughly twice the number of pixels of the full-width half 
maximum of an unresolved line). Clearly, Eq. IA2l can be generalized for resolved, broader features 
by increasing m to span roughly t wice the full-width ha lf maximum of the anticipated broad line. 
The aperture summation method of lLanzetta et al.1 (|1987l ) assigns equal weight to all pixels included 
in the summation around Aj. This can result in an overestimate of the equivalent width detection 
threshold of the data. 



Schneider et al.l (|1993l ) properly treat the relative weighting of pixels adjacent to Aj by weight- 
ing these pixels by the instrumental spread function, ISF. The ISF is then a probability weighting 
function of the flux centered at Aj with pixel weighting, Pj, defined such that 



1, 



(A3) 



where m = 2j Q + 1, and where j Q is a positive integer. For unresolved features, the ISF, and 
therefore the relative values of the Pj, can be taken as a Gaussian function characterized by the 
Gaussian width ctisf = AAj/2.35, where AAj = Xi/R is the full-width at half maximum of the ISF 
centered at Aj for a spectrograph with resolution R. The value of j is chosen appropriately for the 
ISF such that Pi = P m ~ (effectively making sure that the tails of the probability function are 
sampled out to where the probabilities vanish). 

The 1 a equivalent width detection threshold for unresolved features is then computed using 



0£w(Aj) = Di 



Y,P]E 2 {\ k )/I 2 c {\ k ) 

3=1 



- 1/2 


m 











-1 



(A4) 
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where k = i + (j — 1) — j is the "convolution index". The elements of the probability weighting 
function are given by 

Pi = lP~, (A5) 

J = l 

where 

., = «p{-^}. (A6> 

If the feature is redshifted, then Eq, IA4l must be divided by the factor 1 + z to obtain the rest-frame 
equivalent width limit. 

We have generalized Eq. IA4I to (1) account for resolved features of known FWHM, and (2) 
account for regions in which problematic sky-line subtraction in the vicinity of Aj renders the data 
less certain than quantified by the uncertainty spectrum, E(X). We also explicitly include redshift 
dependence. 

To account for resolved features, the <tisf in Eq. IA6I is replaced with the Gaussian width of 
the line spread function (LSF) of the redshifted resolved feature, which is 



0"LSF 



/FWHM(l + z)\ 2 2 
V 2^5 ) + " ISF 



1/2 

(A7) 



where FWHM is the rest-frame full-width at half maximum of the feature, and z is the redshift. 
We have assumed that both the intrinsic line shape and the ISF are well approximated by Gaussian 
functions. 

To account for large residuals from problematic sky-line subtraction, the error in the flux, 
E(X k ), in Eq. IA4I is replaced by E(X k ), which is determined by the quality of the data at pixel k, 

( \r(X k )\ for |r(A*)| > 3tf(A fc ), 
E(X k ) = { (A8) 
[ E(X k ) for |r(A fc )| < 3£(A fe ). 

where |r(Afc)| is the residual of the flux, I(X k ), with respect to the continuum, 

|r(A fc )| = |J(A fe )-J c (A A )|. (A9) 

Using the residual gives a more conservative estimate of the equivalent width detection threshold. 
The residual of a given pixel is used whenever the flux significantly deviates from the continuum. 
In addition to poor sky-line subtraction, sources of this deviation can be large telluric features, 
blending with absorption features, or a problematic continuum fit. Using the following substitution 
for the normalized flux error simplifies the equations for further analysis, 

Y(X k ) = (AlO) 

lc{X k ) 
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Applying these criteria transforms Eq. IA4I from ISchneider et al.l (|1993l ) into the rest-frame 
equivalent width detection threshold calculation used in this work, 



(TEW(\ 



Di 



3=1 



- 1/2 


m 








J =1 



(1+*) 

The rest-frame equivalent width limits presented in Table [3] are quoted at the 3 a level. 



(All) 



A. 2. Uncertainty Assessment 



To quantify the quality of the equivalent width limits, we estimated the uncertainties in the 
ff£lf(Aj). The equivalent width limit in Eq. IA11I explicitly includes the spectrograph resolution, 
R, the absorption line FWHM, the continuum fit, and the central wavelength of the absorption 
line, A/. Assuming that the uncertainties in these quantities can be estimated and are normally 
distributed, the variances in the equivalent width limits are obtained from 



V„ 



d(JEW 

dR 



-i 2 



5R 



+ 



dFWHM 



1 2 



(5FWHM 



+ 



dlc(Xk) 



1 2 



(A12) 



where 5R is the uncertainty in the resolution, 5FWHM is the uncertainty in the FWHM of the line, 
SI C is the uncertainty in the continuum, and ctx^ctew) is the standard deviation in the equivalent 
width limit due to uncertainty in the wavelength center of the line. The terms in Eq. lA12l for which 
explicit indices appear are computed at the center pixel, i, of the feature. 



A. 2.1. Uncertainty in Resolution 



The uncertainty in the resolution, 5R, can be obtained by measuring the full-width half maxi- 
mum of unresolved sky lines and determining the standard deviation. Applying the chain rule, the 
partial derivative of the equivalent width limit with respect to resolution is 

daEW _ da E w dPj Sctlsf 



dR 



where, 



dP j 

d(JLSF 

dR 



dP< da LSF dR 



D; 



l + z) 



c 1 ' 2 



I AG 



2 C 



D 



Pj (Qj - T/S) , 



1 



A. 

cjlsf^ 3 V2-35 



(A13) 

(A14) 
(A15) 
(A16) 
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and where, 



3=1 



rr 3 
LSF 



(A17) 



B = ±2P t T = ±(^it ex JJ^l 

j= l j=1 LSF I- zcr LSF 

c = £j?i*(a*) * = E-pj-%^ 1 - 

j=l j=l ^ ZCJ LSF 

m 

Again, the terms Pj and Qj are evaluated at the center pixel, i, of the feature, j = j a . Note that 
Qj vanishes at the line center. 



A. 2.2. Uncertainty in FWHM 

In the case where the FWHM has a known uncertainty, 5FWHM, the affect on our equiva- 
lent width limits can be calculated. Again, applying the chain rule, the partial derivative of the 
equivalent width limit with respect to FWHM is 



da E w _ da E w dPj <9ct L sf 



dFWHM dPj da LS F <9FWHM ' 
where daEw/dPj and dPj / '<9<tlsf are given by Eqs. IA141 and IAT51 From Eq. IA71 



(A18) 



daLse _ A+z\ 2 FWHM 
5FWHM V 2 -35 J OLSF ' 



A. 2. 3. Uncertainty in Continuum Placement 



We adopt the method of lSembach et al.1 ()199ll ) to calculate the uncertainty in the continuum, 
6I c (Xk), which is obtained directly from the data by taking the rms of the residuals about the 
continuum, a c , from j = 1 to m and multiplying by 0.5. Numerical simulat ions suggest that th is 



provides a conservative estimate of the errors associated with the continuum ( Sembach et al.l[l99ll ) . 
The resulting uncertainty is 

6I c (X k ) = 0.5 a c (X k ), (A20) 
The partial derivative of the equivalent width limit with respect to the continuum is 



1 Di 



do~EW 

dI c {X k ) ~ 2 (1 + z 



(A21) 



-23- 



where, 

G = E^ 



m 

,2 
3 ' 

5=1 



m 



-2P?Y 2 {\ k )I-\\ k ) for £(A fe ) = £(A fc ), 



(A22) 



* = E 



2P 2 I(A fc )/- 2 (A fc ) (1- J(A fc ) I" 1 (A*)) | for £(A A ) = |r(A fc )|. 



The j elements that are selected for the computation of X depend on whether one uses the flux 
error, .E(Afc), or the residual, r(Afc), (see Eq. IA9|) for a given pixel. The continuum contributes a 
large fraction of the error in the equivalent width limit measurements in regions where the residual 
is large, as can be the case with uncorrected telluric features or problematic sky subtraction. In 
cases such as these we select portions around the feature where the continuum estimate is reflected 
more accurately. 



A. 2. 4- Uncertainty in Central Wavelength 

If there is a known uncertainty in the wavelength center of the line, 5A[, its effect on our 
equivalent width limits can be measured directly from the data. The variance of the equivalent 
width limit with respect to the wavelength center of the line is obtained by calculating the rest 
equivalent width limits over the range AAj = ±5Aj(l + z). Eq. IA11I is used for the rest equivalent 
width limit calculations as before. The resulting variance is 

°%(FEW) = ~77 — 7 E ( a EW n ~ (<?EW)) 2 , (A23) 

where Xf and \7 are the upper and lower wavelengths set by (1 + z)(\i ± 5Xi), respectively. N aEW 
is the total number of equivalent width limit calculations, and ((Jew) is the mean equivalent width 
limit, 

1 A ' 

{<TEW) = T7 E ° EW ^ ( A24 ) 

JV "EW 

This technique more accurately reflects the errors due to uncertainties in wave center because 
it takes the actual data into account. If there is a problematic sky subtraction blended with the 
expected position of the feature the calculation will reflect this with a noticeably higher uncertainty. 
Also, if there is an uncertainty in the redshift, z, this can be incorporated into the calculation. 
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Table 1. DLAs in Sample 



DLA 




QSO 


z abs 


Af(Hi)/10 2[) 
[atoms cm -2 ] 


[Zn/H] 


[Fe/H] 


Refs. 


1 


AO 


0235+164° 


0.524 


50+10 






1.2 


2 


Q 


0738+313 


0.091 


15+2 


< -1.14 


-1 63+ 013 


3,4,8,9 


3 






0.221 


7.9+1.4 


< -o.7olg;« 




3,4,8,9 


1 


B2 


0827+243 


0.518 


2.0+0.2 


< +0.30 


-1.02 + 0.05 


4,5,8,9 


5 


PKS 


0952+179 


0.239 


21+3 


< -1.02 




4,5,9 


6 


PKS 


1127-145 


0.313 


51+9 




> -2 


1,6,9 


7 


Q 


1229-020 


0.395 


5.6+1.0 


-0.47 


< -1.32 


7 



Note. — All upper limits are calculated to 3 a. 

" Metallicity for AO 023 5+164 is estimated via X-ray spectroscopy t o be Z = 0.24 + 0.06Zp) ( Chan- 
dra, iTurnshek et aljliopj) and Z = 0.72 ± O.28Z (ASCA & flQg^T. ljunkkarinen et alj|2004l . 

References. — 1. Turnshek et al. (2003); 2. Junkkarinen et al. (2004); 3. Rao & Turnshek (1998); 
4. Kulkarni et al. (2005); 5. Rao & Turnshek (2000); 6. This work; 7. Boisse et al. (1998); 8. Khare 
et al. (2004); 9. Kanekar & Chengalur (2003) 



Table 2. Journal of Observations 









Grating/ 


Slit 




Exposure 


S/N 


Wavelength 




Dispersion 




QSO 


Facility 


Grism 


Width 


Date [UT] 


Time [s] 


[pixel" 1 ] 


Coverage [A] 


Resolution 


[A pixel" 1 ] 


AO 


0235+164 


VLT/FORS2 


GRIS 600z 


1.0" 


2005 Jul 20 
2005 Sep 6 
2005 Oct 1 


8400 


60-150 


7318-10,744 


1880 


1.59 


Q 


0738+313 


APO/DIS 


HIGH 


1.5" 


2004 Dec 15/19 

2005 Feb 4 


40,600 


57-84 


4367-7817 


2040 


0.62/0.84 
(Blue/Red) 


B2 


0827+243 


Kcck/HIRES 




C2/0.861" 


1998 Dec 22 


22,500 


27-114° 


5185-9234 


43,000 


0.04 


PKS 


0952+179 


VLT/FORS2 


GRIS 600RI 


1.0" 


2005 Apr 4/5/7/8 


9000 


57-95 


5298-8622 


1650 


1.63 


PKS 


0952+179 


WHT/ISIS 


R600R 


1.0" 


2004 May 20 


4500 


20-22 


6312-8114 


3790 


0.45 


PKS 


1127-145 


VLT/UVES 


346,564 


1.0" 


2002 Jul 17/18 


24,900 


3S b 


3041-6809 


45,000 


0.05 


PKS 


1127-145 


Gemini/GMOS-S 


R400 


1.0" 


2004 Jun 19 


3600 


38-84 


5962-9998 


960 


2.75 


Q 


1229-020 


VLT/FORS2 


GRIS 600z 


1.0" 


2005 Apr 9/13 


9800 


50-62 


7464-10,000 


1860 


1.59 


Q 


1229-020 


WHT/ISIS 


R600R 


1.0" 


2004 May 20 


3600 


5-13 


7212-9018 


4980 


0.45 



"The low S/N is measured in the red, and the high S/N is measured in the blue. 

b The S/N is measured at the location of the A4428 DIB. The other DIBs are not covered in this spectrum. 



Table 3. Equivalent Width Detections & Limits of DIBs in DLAs 



Rest-Frame EW Detections & Rest-Frame 3 u EW Limits [mA] 



DLA 




QSO 


z abs 


Facility 


A4428 


A5705 A5780 


A5797 


A6284 


A6613 


1 


AO 


0235+164 


0.524 


VLT/FORS2 


741.5+26.2" 


63.2±8.7 6 216+9 6 


<118 bc 


<128 bc 


<95 6c 


2 


Q 


0738+313 


0.091 


APO/DIS 


<115(8) 


<88(6) c 


<74(3) 


<106(3) 


<414(51) 


3 






0.221 


APO/DIS 


<151(11) 


<85(3) 


<78(3) 


<240(16) c 




4 


B2 


0827+243 


0.518 


Keck/HIRES 


<94(14) a 


<32(1) 


<21(1) 






5 


PKS 


0952+179 


0.239 


VLT/FORS2 


<259(19) 


<90(4) 


<289(14) c 


<102(4) 


<108(6) 








0.239 


WHT/ISIS 




<158(4) 


<237(18) 


< 187(3) 




6 


PKS 


1127-145 


0.313 


VLT/UVES 


<68(5) d 
















0.313 


Gemini/CMOS 




<6106(77) cd 


<3162(116) cd 


<341(l) d 


<380(l) d 


7 


Q 


1229-020 


0.395 


VLT/FORS2 




<129(5) 


<115(6) 


<162(5) 


<135(7) 








0.395 


WHT/ISIS 




<292(10) 


<198(8) 


<627(14) c 





"Spectrum provided by Junkkarinen et al. (2004). 
'Spectrum provided by York et al. (2006). 
c Problematic sky subtraction or large atmospheric band. 

^Uncertainty in resolution is unknown and not included in error assessment. 
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Table 4. Model Predictions of DIBs in DLAs 





QSO 


Observed / 




Rest-Frame W Detections, Limits (3 a), & 


: Model Predictions [mA] 


DLA 


ZDLA 


Models 




A4428 


A5780 


A5797 


A6284 


A6613 


1 AO 


0235+164 


Observed 




741.5+26.2" 


216+9 6 


<118 b 


<128 b 


< 9b b 




z = 0.524 


EW[iV(HI)] 






1059 


317 


1572 








EW[B(B - V) = 


0.23] 


482 


104 


43 


368 


50 






EIB - V)r 




U.oO 


n a q 
U.4o 


U.D4 


U.UD 


n a a 
U.44 


2 Q 


0738+313 


Observed 




<115(8) 


<88(6) 


<74(3) 


<106(3) 


<414(51) 




z = 0.091 


EW[A r (HI)] 






247 


89 


532 








EW[E(B - V) = 


0.04] 


84 


18 


8 


91 


3 






EiB - V)tirr, 




0.05 


0.20 


0.40 


0.05 


1.91 


3 Q 


0738+313 


Observed 




<151(11) 


<85(3) 


<78(3) 


<240(16) 






z = 0.221 


EW[A r (HI)] 






114 


46 


299 








EW[B(B -V) = 


0.04] 


84 


18 


8 


91 








EiB - W 




0.07 


0.19 


0.42 


0. 14 




4 B2 


0827+243 


Observed 




<94(14) 


<32(1) 


<21(1) 








z = 0.518 


EW[A r (HI)] 






22 


11 










EW[B(B - V) = 


0.04] 


84 


18 


8 










EIB - 




0.04 


0.07 


0.11 






5 PKS 


0952+179 


Observed 




<259(19) 


<90(4) 


<237(18) 


<102(4) 


<108(6) 




z = 0.239 


EW[AT(HI)] 






370 


127 


720 








EW[B(B - V) = 


0.04] 


84 


18 


8 


91 


3 






E{B - V) lim 




0.12 


0.20 


1.29 


0.05 


0.50 


6 PKS 


1127-145 


Observed 




<68(5) 


<6106(77) 


<3162(116) 


<341(1) 


<380(1) 




z = 0.313 


EW[AT(HI)] 






1086 


324 


1601 








EW[B(B - V) = 


0.04] 


84 


18 


8 


91 


3 






E(B - V) lim 




0.03 


14.13 


17.73 


0.21 


1.75 


7 Q 


1229-020 


Observed 






<129(5) 


<115(6) 


<162(5) 


<135(7) 




z = 0.395 


EW[AT(HI)] 






75 


32 


219 








EW[B(B -V) = 


0.04] 




IS 


8 


91 


3 






E(B - V) nm 






0.29 


0.62 


0.08 


0.62 



"Spectrum provided by Junkkarinen et al. (2004). 
^Spectrum provided by York et al. (2006). 

c The limits are those that are the best constrained from Table [3] The model EW[A r (HI)] refers to the predicted equivalent 
widths (mA) from the Galactic best— fit lines in Figure [9] The values used for the N (HI) are taken from Table [T] The model 
EW[£(B - V) = 0.23] re fers to the predicted equ ivalent widths (mA) from the Galactic DIB-E(B - V) best-fit lines with the 
known E(B - V) = 0.23 jjunkkarinen et aljISoM) . The model EW[B(B - V) = 0.04] refers to the predicted equivalent widths 
(mA) from the Gal actic DIB— E(B — V) best-fit lines with the upper limit of E(B — V) = 0.04 as measured in high-z DLAs by 
lEUison et al] ||2005). Fo r all systems, the E(B - V) = 0.04 except for AO 0235+164 which has a measured E(B - V) = 0.23 
Tjunk1tarrn^ri~et**aL 120041) . The value E{B — V)n m is the E{B — V) upper limit inferred from our equivalent width limits and the 



Galactic DIB-£(S - V) best-fit lines for the A5780, A5797, and A6284 DIBs l lWeltv et al.ll2006ri and the A4428 and A6613 DIBs 
(unpublished, T.P. Snow, private communication). 
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Table 5. Adopted Reddening and Gas-to-Dust Ratios 



DLA 




QSO 


%abs 


E(B - V) 
[mag] 


gas/dust 
[cm -2 mag -1 ] 


Contraint 
DIB 


1 


AO 


0235+164 


0.524 a 


0.23 


19.2X10 21 




2 


Q 


0738+313 


0.091 


< 0.05 


> 30 xlO 21 


A6284 


3 






0.221 


< 0.14 


> 5.6X10 21 


A6284 


4 


B2 


0827+243 


0.518 


< 0.07 


> 2.9xl0 21 


A5780 


5 


PKS 


0952+179 


0.239 


< 0.05 


> 42 xlO 21 


A6284 


6 


PKS 


1127-145 


0.313 


< 0.21 


> 24 xlO 21 


A6284 


7 


Q 


1229-020 


0.395 


< 0.08 


> 7.0xl0 21 


A6284 



"Values quoted from Junkkarincn ct al. (2004). 
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Fig. 1 — Equivalent width limit analysis of the Q0738+313 z = 0.091 DLA. —(a) Relative flux in 
the region of the A5780 and A5797 DIBs (vertical long dashed line and vertical long dashed-dotted 
line, respectively). The histogram is the relative flux data, and the short dashed lines are the 
associated ±3 a flux errors. — (6) Relative sky flux. — (c) The 3 a rest equivalent width limits 
in mA. The two peaks in panel (c) are due to the problematic sky subtraction from the [OI] sky 
emission line at A ~ 6300 A and the poor continuum fit at A ~ 6314 A where the residuals are 
used, as explained in Eq. IA81 
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Fig. 2.— VLT/F0RS2 spectrum of the z = 0.524 DLA toward AO 0235+164. Plotted are the 
expected positions of the (a) A5780/5797, (b) A6284, and (c) A6613 DIBs. The upper sub-panels are 
the normalized flux of the data and models. The center sub-panels are the sigmas of the associated 
data fluxes normalized by the continuum. The lower sub-panels are the sky fluxes normalized by the 
continuum. The smooth curves are model predictions (see text). The thin curves are the expected 
DIB profiles given the measured DLA iV(Hi) and the Galactic DIB-iV(Hi) r elation. The thick 



curves are the expected DIB profiles using the measured E(B — V)= 0.23 from I Junkkarinen et al 
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Fig. 3.— Same as Fig. El but for the APO/DIS spectrum of the z = 0.091 DLA toward Q0738+313. 
Plotted are the expected positions of the (a) A4428, (b) A5780/5797, (c) A6284, and (d) A6613 DIBs. 
The smooth curves are model predictions (see text). The thin curves are the expected DIB profiles 
given the measured DLA iV(Hi) and the Galactic DIB-iV(Hi) relation. The thick curves for this 
DLA are the expec ted DIB profiles using the E{B — V) = 0.04 upper limit for high-z DLAs from 



Ellison et al.1 fl200J ). 
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Fig. 4.— Same as Fig.H but for the APO/DIS spectrum of the z = 0.221 DLA toward Q0738+313. 
Plotted are the expected positions of the (a) A4428, (6) A5780/5797, and (c) A6284 DIBs. 
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Fig. 5. — Same as Fig. [3l but for the Keck/HIRES spectrum of the z = 0.518 DLA toward 
B2 0827+243. Plotted are the expected positions of the (a) A4428 and the (6) A5780/5797 DIBs. 
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Fig. 6.— Same as Fig. H but for the VLT/FORS2 and WHT/ISIS spectra of the z = 0.239 DLA 
toward PKS 0952+179. Plotted are the expected positions of the (a) A4428, (b) A5780/5797, (c) 
A5780/5797 (WHT/ISIS), (d) A6284, and (e) A6613 DIBs. The WHT/ISIS spectrum provides the 
adopted limit for the A5797 DIB; the VLT/FORS2 spectrum provides the adopted limits for the 
remaining DIBs. 
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Fig. 7. — Same as Fig. [3J but for the Gemini/GMOS-S spectrum of the z = 0.313 DLA toward 
PKS 1127-145. Plotted are the expected positions of the (a) A5780/5797, (b) A6284, and (c) A6613 
DIBs. 
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Fig. 8.— Same as Fig. O but for the VLT/FORS2 spectrum of the z = 0.395 DLA toward Q1229- 
020. Plotted are the expected positions of the (a) A5780/5797, (6) A6284, and (c) A6613 DIBs. 
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Fig. 9.— The DIB equivalent width-iV(HI) relations (jWeltv et al.1 12009 ) with our DLAs added, 
-(a) A5780 DIB. —(b) A5797 DIB. —(c) A6284 DIB. The solid lines are the best-fit weighted 
Galactic lines. The region enclosed by the dotted lines contains the Galactic data. The region 
enclosed by the dashed lines contains the LMC data. The region enclosed by the dot-dash lines 
contains the SMC data. Error bars are 1 a, and upper limits are marked with arrows. The vertical 
error bar s for AO 02 3 5+16 4 in panel (a) are smaller than the point size and all values for this DLA 
are from lYork et al.l ( 2006 ). 
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Fig. 10. — The logarithm ratio of the measured equivalent width limits to the equivalent width 
predicted for multiple values of reddening for each DLA. — (a) A4428 DIB. — (b) A5780 DIB. — (c) 
A5797 DIB. —(d) A6284 DIB. —(e) A6613 DIB. The equivalent widths of the A5780 A5797, and 



A6284 DIBs, for each reddening, were computed using the best-fit Galactic lines from lWeltv et al 



(|2006l ). The equivalent widths of the A4428 and A6613 DIBs, for each reddening, were computed 
from the relations in Eqs. [J] and [5] (unpublished, T.P. Snow, private communication). The DLAs 
labeled are as follows: (1) AO 0235+164, (2) Q0738+313, z = 0.091, (3) Q0738+313, z = 0.221, 
(4) B2 0827+243, (5) PKS 0952+179, (6) PKS 1127-145, (7) Q1229-020. We leave the panels of 



the A4428 and A5780 DIBs for AO 0235+164 blank because they ar e detections flj unkkari nen et al 
20041 : lYork et alJbood ) with known reddening, E(B - V) = 0.23 (jjunkkarinen et al.ll20o3 l. The 



zero point (horizontal dotted line) on each plot marks the estimated maximum reddening of the 
DLA determined from the measured equivalent width limit, or E(B — V r )n m (see Table H]). The 
vertical error bars are the 1 a errors for the E(B — V)ii m based on the equivalent width limit 

nnpprtnintip« and thp lmpprtnintips in trip FlTR-rpdHpnincr law «lnnps in T?<n<i fTl [9l Rl [Zl and [Kl Thp 
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Fig. 11. — The adopted upper limits on the gas-to-dust ratios of the DLAs in our sample. For 
comparison, the various lines provide ranges of the measured values (given in the text) for the Milky 
Way (MW), Large Magellanic Cloud (LMC), and Small Magellanic Cloud (SMC). The limits are 
illustrated with solid leftward arrows, based upon the A5780 and A62 84 DIBs. The AO 0235+164 
reddening measurement (star datum point) is from I Junkkarinen et al.l (|2004l ). 



